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Abstract 

We define the generalized-Euler-constant function 7(2) = XmLi z n_1 (- — log 1 ^) 
when \z\ < 1. Its values include both Euler's constant 7 = 7(1) and the "alternating 
Euler constant" log - = j(— 1). We extend Euler's two zeta-function series for 7 to 
polylogarithm series for 7(2). Integrals for 7(2) provide its analytic continuation to C — 
[l,oo). We prove several other formulas for 7(2), including two functional equations; 
one is an inversion relation between 7(2) and 7(1/2). We generalize Somos's quadratic 
recurrence constant and sequence to cubic and other degrees, give asymptotic estimates, 
and show relations to 7(2) and to an infinite nested radical due to Ramanujan. We 
calculate j(z) and "y'(z) at roots of unity; in particular, 7 ; (— 1) involves the Glaisher- 
Kinkelin constant A. Several related series, infinite products, and double integrals are 
evaluated. The methods used involve the Kinkelin-Bendersky hyperfactorial K function, 
the Weierstrass products for the gamma and Barnes G functions, and Jonquiere's relation 
for the polylogarithm. 
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1. Introduction 

In this paper we introduce and study the generalized- Euler- const ant function j(z), defined 

by 

1 ( z \ = z n-l (}_ _\ og U+1 \ = f f 1 ~ X dxdy ^ 

^ \ n n J Jo Jo {1- xyz){- log xy) 

where the series converges when \z\ < 1, and the integral gives the analytic continuation 
for z G C — [l,oo). The function 7(2) generalizes both Euler's constant 7 = 7(1) and the 
"alternating Euler constant" log - = 7(— 1) [T%] . where 7 is defined by the limit 

7 = lim ( 1 + - + ••• + -- logn J = 0.57721566 .... (1) 

n-»oo y 2 n J 

In Section 2, we extend Euler's two zeta-function series for 7 to polylogarithm series for 7(2) 
(Theorem Q); we also show another way in which the function 7(2) is related to the extended 
polylogarithm (Theorem 0J). We give a variant of the series definition of 7(2) which is valid on 
a different domain (Theorem |SJ). We prove two functional equations for 7(2) (Theorems |H1 and 
|7J); the first is an inversion formula relating 7(2) and 7(1/2:), and is proved using Jonquiere's 
relation for the polylogarithm. 

In Section 3, we generalize one of Somos's constants, and show relations with the function 
7(2) (Theorem^ and with an infinite nested radical due to Ramanujan (Corollary |2J). We also 
generalize Somos's quadratic recurrence sequence to cubic and other degrees (Theorem EJ) and 
provide asymptotic estimates (Lemma ^ and Theorem I10J). 

In Section 4, we calculate the value of 7(2;) at any root of unity ( Theorem [T3J); the proof uses 
a result based on the Weierstrass product for the gamma function ( Theorem IT2*j). Using the 
Kinkelin-Bendersky hyperfactorial K function, in Section 5 we compute the derivative ^'(00) at 
a root of unity u 7^ 1 (Theorem EJ); the proof involves a summation formula derived using the 
Barnes G function (Theorem Ho]) . In particular, we show that 7'(— 1) involves the Glaisher- 
Kinkelin constant A (Corollary 0J), and that it is related to an infinite product (Example fTTj) 
essentially due to Borwein and Dykshoorn 0. 

Other infinite products occur in Corollary Q Example 0] (where we generalize an accelerated 
product for pi 17J), Section 3, and Example We evaluate some new double integrals in 
Examples IH1 and ITUl Theorem IT41 and Corollary 0] 
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2. The Generalized-Euler-Constant Function 7(2) 



Denote by Z + , R, and C the sets of positive integers, real numbers, and complex numbers, 
respectively. If z £ C and 2 / 0, define 

logz = In |z| + i Argz (— 7r < Arg z < 71). 

The limit definition of Euler's constant is equivalent to the series formula 



7 = Yl ( - - iog^^ ) 

z — ' V n n J 

n=l v 7 



(2) 



(see [16J). The corresponding alternating series gives the "alternating Euler constant" [IS], JH] 
(see also [Zj) 

log 1 = y^(-l)"- 1 ( - - bg = 0.24156447 .... (3) 

7r V n n / 

n=l v 7 

The main subject of this paper is the following function j(z), which generalizes ^ and (0). 

Definition 1 The generalized- Euler- constant function j(z) is defined when jz] < 1 by the 
power series 



= (I- log i±±) 

n=l x ' 



(4) 



which converges by comparison to series ^ for 7. 



Example 1 In addition to 7(1) = 7 and j(— 1) = log-, Definition ^ gives 7(0) = 1 — log 2. At 
2 = 1/2, the function takes the value 



7 (-)=21og-, (5) 



where 



a = y 1 V 2v / 3^ = lV22 1 /4 3 i/8 ... = JJ n i/2" = 1.66168794 . . . 

n=l 

is one of Somos's quadratic recurrence constants [Uj, ^3] Sequence A112302] (see also 
and [3 p. 446], where the notation 7 is used instead of a). To see this, write (J3J) with z — 1/2 

as 

/1\ / 1 21og(n+l) logra 

Tin " 



2 ^ \ 2 n n 

n=l 



2 

2(log2-21oga + log(x) =21og-. 

a 
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Euler gave two zeta-function series for his constant 7 (see, for example, equations 3.4(23) 
and 3.4(151)]), 



C(fc) - 1 

k ' 



(6) 



k=2 k=2 

We generalize them to polylogarithm series for the function 7(2). 



Theorem 1 If \z\ < 1 and z/ Li^(^) denotes the polylogarithm 4., Section 1.11], [7], [T7j . 

defined for k = 2, 3, . . . fry £/te convergent series 



Li ^) = E^ 



n=l 



then 



(7) 



k=2 



\A k (z) - z 



If in addition z 7^ 1, then 

00 

z 2 7 (z) = z + (1 - z) log(l -z)-J2 

k=2 

Proof. If \z\ < 1, then Definition ^ and the expansion 

io g (i-™) = -f]^ (H<i,^^i) 

give 

z 1 (z) = z(l-log2) + J2 



m=l 



(9) 



n=2 







- - log 


K)] 


n 





00 00 



fc=2 n=2 fc=2 



/cn fc 



It is easy to see that the double series converges absolutely, so we may reverse the order of 
summation, obtaining 



fc=2 



n=2 



fc=2 



A' 



This proves (J7J). 

Now take z ^ 1 with |z| < 1. Definition Q and formula © imply that 



00 „ 00 

z 



7 (z) = * E - - E ] °g — = -* lo g(! - *) - E lo S ^7 



n=l 



n=l 



n=2 



where we have re-indexed the last series. Thus, using © again, we can write 

oo „ oo 

z 2 7 (z) = -zlog(l - z) +log(l - z) + + ^z n lo, 



n * — ' n 

n=l n=2 



00 n f i\ 
;i- z )io g (i-z) + z + V> - + io g 1-- 

z — ' n \ n J 



n=2 

The last series is equal to the absolutely convergent double series 

\ ^ \ ^ Z \ ^ l \ ^ Z _ \ ^ ^k{z) — Z 

^ h n k J- 2-^i n k k 

n=2 k=2 k=2 n=2 k=2 

and the proof of (JHJ is complete. □ 

Example 2 Let z = 1 in (JJJ), and let z tend to 1~ in (jHJ). Since 7(1) = 7 and Li^(l) = C(fc), 
we recover the two zeta series (JBJ) for Euler's constant. 

Now take z = -1, and substitute 7(-l) = log^ and Li fc (-1) = (2 1_fe - l)C(fc) (see [6, 
Section 9.522], [Zj, [T2j, [Tj3 Chapter 3]). Formula (J2J) and the first equality in (jHJ) give the zeta 
series (compare [TBI equation (4)] and [213 equation 3.4(25)]) 



k=2 ' k=2 



Simplification in formula (jHJ) yields a zeta series for the logarithm of pi: 

l-(l-2 1 - fc )C(fc) 



log 7T = 1 + ^ 



k=2 



k 



Finally, set z = 1/2. Using formula (JEJ) for the Somos constant a, and defining Li^z) 
— log(l — z), we can write formulas (J7J) and (jHJ as 



fc(l/2)-l 



fc=i 

Theorem 2 T/ie function j(z) is continuous on the closed unit disk 

D = {z e C : |z| < 1}, 

and holomorphic on the interior of D. However, the left-hand derivative at z — 1 does not 
exist; more precisely, 

lim 7(1) ' 7(t) = +oo. (10) 
t-*i- 1 — £ 

In particular, the Taylor series expansion ofj(z) at z = /ias radius of convergence 1. 
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Figure 1: The image of the unit circle under the function 7(2;) 



Proof. For z G D, series (JIJ is majorized by series (J2J) for Euler's constant. Therefore, series (jlj) 
converges to 7(2) uniformly on D. It follows that 7(2) is continuous on D, and holomorphic 
on the interior of D. 
If < t < 1, then 



7(1) "7(0 = D 1 - ^ 



n=l 

00 

Ed 

n=l 



1 , / 1 

- - log 1 + - 

n \ n 



t 



t -n (±_ 



\ 2n 2 3n 3 



n 3 4n 4 / ~ ^ V2n 2 3n 3 / 

7 n=l x ' 



for N = 1, 2, 3, . . .. Multiplying by (1 — £) 1 , and letting t tend to 1 , we obtain 

N / , ^ \ TV 



lim inf 

t-»l- 



7(1) ~l(t) 
l-t 



n=l v 7 n=l v 



n 6n 2 



The last sum tends to infinity with N, and the theorem follows. 



□ 



The image of the unit circle under the function j(z) is shown in Figure 
In [13 HEl HZ| , the first author represented Euler's constant 7 = 7(+l) and the alternating 
Euler constant log - = j(— 1) by the double integrals 

-1 /■! 



7 (±1) 



1 — x 



■ dx dy. 



'0 J i.^Txy)(-\ogxy) 

We extend this to integrals for the function 7(2), and obtain its analytic continuation to the 
domain C — [1, 00). 
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Theorem 3 If \z\ < 1, then 



1 rl 



o Jo 



x 



l{ z ) = / / T\ \7 — i rdxdy= / — — dx. (11 



'1 — xyz)(— \ogxy) J (1 — xz)\ogx 



1 1 — x + log x 



The integrals converge for all z G C — (l,oo), and provide the analytic continuation of the 
generalized- Euler- constant function j(z) for z G C — [1, oo). 

Proof. By |7J Theorem 4.1], if z G C — (1, oo) with z ^ 0, and if 3?(s) > —2 with s ^ —1, then 
for u > 



(1 — z)$(z, s + Lit) - u 



Here $ is the Lerch transcendent [41 Section 1.11], Section 2], the analytic continuation of 
the series 

d>(z, g ,n) = ]P (13) 

n=0 ' 

which converges for all complex s and all u > when |z| < 1. In particular, $(z, 0,u) = 
1 + z + z 2 -\ = (1 - z)- 1 . It follows, letting s tend to -1 in (fP2j) . that 

Jo Jo (l-arj/z)(-logxj/) z ds z 

Now set it = 1 and let I(z) denote the double integral in (JTJJ. Using (|T3*Jl . and re- indexing, we 
obtain, when < |z| < 1, 

(15) 

OO 

E - — E * n_1 lo §^ + x ) + E * n_1 lo s n = 7(z)> 



n=l n=l 



n 

n=l n=l n=l 



by Definitional 

This proves that 7(2;) = J(z) when < \z\ < 1. It follows, since the functions 7(2) and J(z) 
are continuous on D — {1} (in fact, on D), that 7(2) = I(z) also when |z| = 1 z and when 
z = 0. Finally, 7(1) = 7 = 1(1), from [E]. 

It remains to show that I(z) is equal to the single integral in To see this, make the 

change of variables x = X/Y, y = Y, and integrate with respect to Y. □ 

The graph of the function 7(2) for z G (—00, 1] is shown in Figure El Its properties are 
easily verified from (fTTj) and (|10jl . Namely, /or 21 reaZ, £/ie graph 0/7(2;) positive, increasing, 
and concave upward; it is asymptotic to the negative real axis, that is, lim^oo 7(— z) = 0; and 
the tangent line at the point (1,7) is vertical. 
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—14 —12 — 1 —8 —6 —4 —2 

Z 



Figure 2: The function 7(2) for real z 



Example 3 Using Theorem El and formula (0) for Somos's constant a, we recover the evalua- 
tions from [7j 



Jo 



dxdy = / — — cfe = log a. 



(2- xy)(-\ogxy) J (2-x)(-\ogx) 

To see this, set z = 1/2 in (fTTj) . and replace 7(1/2) with 21og(2/cr). Now multiply by 1/2, and 
subtract the result from the equations 



1,1 x ,1 j 



'0 jo 

which are easily verified 



(2-a^/)(-logxj/) Jo 2-x 



dxdy = / c?x = log 2, 



The function 7(2) is also related to the polylogarithm in a different way from that in 
Theorem ^ 

Theorem 4 For all z G C — [1, 00) the relation 

z 2 7 (z) = (1 - z)Li (z) - z log(l - z) (16) 

holds. Here the prime ' denotes d/ds, and Li s (z) is the extended polylogarithm, the analytic 
continuation |U Section 1.11], Section 2], [T7| Section 5] of the series 

00 

z 



L Uz) = £-- (17) 



1 n 
n=i 



If \z\ < 1, then the series converges for any seC, and the relation can be written 

oo 

z-f(z) = - log(l -z)-(l-z)J2 ^ lo S n - ( 18 ) 



n=l 



We give two proofs. 

Proof 1. Multiply the last two equations in ()15|) by z, with \z\ < 1. Using ©, formula (fTBj) 
follows. Now multiply (fTHj) by z. Using (fTTj) . we obtain relation ()16|) . and the theorem follows 
by analytic continuation. □ 

Proof 2. Setting u — 1 in (fPfj) . we may replace the integral with 7(2), by Theorem El Now 
multiply the equation by z 2 . Using the formulas — log(l — z) — z§(z, 1, 1) (from Q and (fTTHO 
and 

Li.(z)=z$(z,s,l) (19) 

(from (HH) and (HSI)), relation (JH3) follows. If \z\ < 1, substituting IflTjl in (fTTH) yields (HHJ). □ 

An application of Theorem 0] is a formula for 7(2:) involving a series which converges if 
U(z) < 1/2, when |^| < 1. 



00 , \ 11 n , \ 

mi - + E E(-!) fe+1 U m* + 1). 

n=l ^ ' k=0 ^ ' 



Theorem 5 < 1/2, £/*en 

zj(z) = - io g (i - z) + > ; ( — ) > ;(-i) fe+i ( ) io g (* + 1). (20) 

Equivalently, if \w\ < 1, then 

^ {t~) = - w ) ~ «0 + (1 " «0 £ W " E^ 1 )" (?) lQ g( fc + ^ ( 21 ) 

^ 1,7 ^ n=l fc=0 ^ ' 

Proof. By a special case of Corollary 5.1], if 9ft(z) < 1/2, then the double sum in (|2U|) is 
equal to the product (1 — z)^{z, 0, 1). Using (JT§|| and ((TOJ) , the first statement follows. To 
prove the second, note that if |iu| < 1 and z = f^-, then < 1/2, and PU]) implies (}2*T]) . □ 

Corollary 1 For m = 1, 2, 3, . . . , define the infinite product 

oo / n 

4=n iK fc+i ) H)W0 

rt=l \fc=0 , 

T/ien i/ie product converges, and 

7 (-m) = — log— — . 
m P m 
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Proof. When z = —m, the series in (j2"Uj) is equal to logP m . □ 

Example 4 Let m — 1. The product Pi is the acceleration of Wallis's product for tt/2 in J7| 
(see also (Zj): 



Pi 

Thus 



2 \l/2 / 2 2 \V4 / 2 3.4\ 1 /8 / 2 4.44 \V16 ^ 



1/ V 1 - 3 / V 1 ■ 3 3 / V 1 ■ 3 6 ■ 5 



2 4 
7 (-l) = log — = log-, 



confirming the value of 7(— 1) in Example ^ 
With m = 2, we get 

7(-2) = -log-, 

where 

- 2 \ 2/3 / 2 2 \ 4/9 4 \8/27 / ^ ^ N 16/81 

P 2 = ( - ] ' - ^ ' " ' ' 



i; vi ■ 3/ vi-3 3 ; u-3 6 -5 

We now prove a functional equation for the generalized- Euler-const ant function which ex- 
presses 7(1/2) in terms of 7(2). 

Theorem 6 The following inversion formula holds for all z G C — [0, 00) with ^s(z) > 0: 



7 ( - ) = z 3r y(z) — Triz + z 2 log(l — z) — zlog h z(l — z) 

\zj z 

where ip(z) is the digamma function 

d[\ogT(z)] _ T'(z) 



7li /log .2 

7 + log 27T + — + ip 



2 T \ 2iri 



dz r(z) 



Proof. The extended polylogarithm Li s (;z) (see Theorem 0} satisfies Jonquiere's relation jU 
Section 1.11], 

(2i) s e™/ 2 / log z 



U s {z) + e™U s {l/z) = v r(g) C (^1 - a, ^-J , (22) 

where ((s,w) is the Hurwitz (or generalized) zeta function |H Section 1.10], the analytic con- 
tinuation of the series 



((s,w) = V7 — : n 



n=0 
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which converges if dt(s) > 1 and w G C — {0, —1, —2, . . .}. If \z\ < 1 and s — 0, the sum of 
the series (II 7|) for Li s (z) is Li (z) = z(l — z) -1 , which by analytic continuation holds for all 
complex z 1. It follows that at s = the left side of (}22|) is equal to —1. Therefore, the 
derivative of (|2*2"|) with respect to s at s = is 



Li' (z)+Li' (l/z) + 



Til 



Z - 1 



lim - 

s^O S 



(^'e™/ 3 A _ ^ lug : 



r s 



2m 



To compute the limit, we use the Taylor series for 



(2Tr) s e 



s iris/2 



exp 



log 2tt 



Til 



together with the estimates [T^J Equation 43:6:1] 

1 



T(s) 



and H Section 1.11] 



({l-s,w) = if>(w) + 0(s) (3?H>0), 



which are valid for s tending to 0. If Q(z) > 0, the result is 



W (z)+U' (l/z) + 



Til 



Z - 1 



-7 — log 2tx — 



Til 



1> 



logz 
2m 



For z G C — [0,oo), we may use (fTB^l to replace Li^) and Li' Q (l/z) with expressions for 
them involving 7(2) and 7(1/2), respectively. Solving for 7(1/2;), we arrive at the inversion 
formula. This proves it when ^s(z) > 0. Since each term in the formula is continuous on the 
set {z G C| SR(s) < 0, $$(z) > 0}, the theorem follows. □ 

Example 5 Take z — — 1. Using the values log(— 1) = ni and if) (1/2) = —7 — log 4 (from 
Gauss's formula (201 equation 1.2(47)] 



"7 



— cot log q + cos log I 2 sin — 1 



k=l 



. kTi\ 



Q J 



(23) 



where < j < q), we again obtain 7(— 1) = log -. For a related application of Theorem see 
the proof of Corollary 

Remark 1 Setting z = 1/w in Theorem |U1 we obtain an inversion formula valid for w G 
C - [0, 00) with 9f(tw) < 0. 
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The next theorem gives a second functional equation for the function y(z). The equation 
relates the quantities y(z),y(—z), and y(z 2 ). 

Theorem 7 The following reflection formula holds for all z G C — ((— oo, —1] U [1, oo)): 

z(l + z)y(z) + z{\ - z)y{-z) = 2z 3 y{z 2 ) - 2z log 2 + (1 + z) log(l + z) - (1 - z) log(l - z). 
Proof. If \z\ < 1, then using (fTTj) we see that 



1 00 2 l o 

- (Li' (.) + W (-z)) = - J> 2n log 2n = + Li' ( 

n=l 



z 2 ) 



where the prime ' denotes d/ds. The relation between y(z) and Lig(,z) in TheoremEJthen yields 
the desired formula, and the result follows by analytic continuation. □ 

Remark 2 Theorem [7| can be generalized, as follows. Given an integer q > 1, we choose a qth 
root of unity 1, and obtain the average 

- 2^ Ll o(^) = -^—7 + Llo(2 

y j=0 

Theorem |3] then translates this into a formula relating y(z),y(uz), . . . , 7(o;' ? ~ 1 z) and y(z q ), valid 
for all 2 G C such that u^z G" [1, oo) for j = 0, 1, . . . , q — 1. 



3. A Generalization of Somos's Quadratic Recurrence Constant 

We begin this section by generalizing both Somos's quadratic recurrence constant a and 
its relation (jJJ) with the function y(z). By convention, if a > and t > 0, we assume that 
= a 1 /* > 0. 

Definition 2 For t > 1, the generalized Somos constant o~ t is given by 

n=l 

The convergence of the infinite product for 07 follows from the convergence of the series 

^logn _., /l\ 19$ (\ n \ 

n=l x ' x ' 

where $ is the Lerch transcendent (|13p. Note that for £ = 2 we get Somos's constant <r = o~ 2 - 

The following result relates the generalized Somos constant at to the function y(z), essen- 
tially generalizing Example ^ 
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Theorem 8 Fort > 1, the generalized Euler constant 7 '(1 ft) and the generalized Somos con- 
stant a t satisfy the relation 

*log - I (25) 



7 



it - l)o* t 



In particular, 



lim taj , 1 = e 



(26) 



Proof. To prove ()25jl . set 2 = 1/t in (|18j) and use (|24jl . (The case t = 2 is proved in Example [I]) 
To prove 1)26)1 . replace £ with £ + 1 in ()25j) . let £ tend to + , and use linx,,^!- 7(2;) = 7(1) = 7. □ 

Remark 3 From Theorem |S] and the reflection formula in Theorem [7[ one can also express 
7(— 1/t) in terms of at and a t z. 

The next theorem generalizes a result due to Somos [H] (see p. 446] and |2*T|). 

Theorem 9 Fzx £ > 1, and define the sequence (g n ,t)n>o via the generalized Somos recurrence 

9o,t = 1, 9n,t = ng l n _ 1:t (n > 1). (27) 

JTiere we /icwe £/ie explicit solution 

"10$ /l 

7 as 

Proof. This follows by induction on n,. □ 

For £ = 2, equation (|77jl is Somos's quadratic recurrence g n = ng\_ x . Somos [T4] . [TBI 
Sequence Al 16603] (see also (3] and j5TJ) gave the following asymptotic formula for g n = g n ^ 
as n tends to infinity: 

g n ~a 2n (n + 2- n" 1 + An' 2 - 21n~ 3 + 138n" 4 - 1091n~ 5 + . . (28) 

We extend this to an asymptotic formula for g n t given any fixed t > 1. 

Lemma 1 For t > 1, let (g n ,t)n>o be the sequence defined by the recurrence ()27j). For x G 
[0,oo), define ft{x) G [1, 00) by the infinite product 



g„ t = a t exp 



0,n + l 



m=l 



,i/t" 



m=l 



TTien /or eac/i N > 2 and n > 1, there exists a positive number fi = fi(N,n,t) such that 



9n,t = afn-^ 



N-l 



fc=l 



1 d k f t 
n k k\ dx k 



(0) + 



1 d N f t 
n N N\ dx N 



(M) 



(29) 



Moreover, for fixed N > 2 and /ixecZ £ > 1, 



lim fi(N,n,t) = 0, lim 
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Proof. Since t > 1 and x > 0, the product for converges. It follows from Theorem El and 
the identity YZ=i r ™ = V(* ~ X ) that 



ft* = ^r^- i/(t - i} /f ' 1 



r- 1 



Thus, by Taylor's theorem with remainder, it suffices to show that ft{x) is infinitely differen- 
tiable with respect to x for x G [0, oo), the derivatives at x = being right-sided. That in turn 
follows (see, for example, [TTl p. 342, Theorem 4]) by noting that for k > 1 the fcth termwise 
derivative of the series for \ogf t (x), namely, 



-l) fc - 1 (A;-l)!y — — 
4 t m (l 



k 

m 



+ mx) k 



m=l 

is uniformly convergent on [0, oo), as it is majorized by (k — 1)! Ylm=i mk t~ m - □ 

Theorem 10 For fixed t > 1, the sequence (g n ,t)n>o satisfies the following asymptotic condition 
as n tends to infinity: 



„ t _ t™ -i/(t-i) 
yn,t — Of n 



t t(t 2 -t-1) t(2t 4 + 1 3 - lit 2 + 7t + 2) n (}_ 

+ (1 - t) 2 n ~ 2(1 -i) 4 n 2 + 6(l-t) 6 n3 + ° \tf 



Proof. By Lemma CJ we only have to compute |^(0). For k = 1, 2, . . . define 0a,. (t) by 
**(') = ^(0) = <-!)"<* " D' E £ = - !) !Li - (l) ■ 

m=l ^ ' 

Using the identity 

& (t) = (fc-l)t^I(t) (fc>l), 

we compute that 

t(l +t) , 2t(t 2 + 4^ + 1) 



1-t) 2 ' rzw (1-t) 3 w (l-t) z 



Since / t = e log ^ and /t(0) = 1, Faa di Bruno's formula for the fcth derivative of the 
composition of two functions (see, for example, [SJ p. 21]) yields 

dx k ' mi! ■ • ■ mfc! \ j! 

where the sum is over all fc-tuples (mi, . . . , mk) of non- negative integers such that 

mi + 2m 2 + 3m 3 + . . . + km,k = k. 
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In particular, 

|^(o) = <Mt), ^( O ) = 1 (t) 2 + 2 (t), ^Ii(o) = M t) 3 + 3Mt)Mt) + Mt)- 

Applying ()29|) with N = 4, the theorem follows. □ 

Example 6 Taking t = 2 gives the first four terms of the asymptotic formula (J28)) for Somos's 
quadratic recurrence sequence g n = g n ^ = ng\_ x 2 = 1,2, 12, 576, 1658880, . . . [T3J Sequence 
A052129]. With t = 3 we get an asymptotic formula [3 Sequences A123853 and A123854] 
for the cubic recurrence sequence g n $ = ng^_ 13 = 1,2,24,55296,845378412871680,... [T3J 
Sequence A123851], namely, 



V 1 — 1 II 
9n,3 = 0" 3 71 



' _3 15_ 113 /J_ 

+ An 32n 2 + 128n 3 + In 4 



In — > oo 



where u 3 is the cubic recurrence constant Sequence Al 23852] 



(T 3 = y iy 2^3^ = lV3 2 V93i/27 . . . = 1.15636268 .... 

The generalized Somos constants at are connected with some of Ramanujan's infinite nested 
radicals and with the Vijayaraghavan-Herschfeld convergence criterion for them, as follows. Fix 
t > 1. Given (t) > for n = 1, 2, . . . , define b n = b n (t) by 



\ 



a x + y a 2 + y a 3 H h y «n-i + v^- (30) 

If the sequence (6 n )n>i converges, we call the limit = b^t), and we define 



a\ + V a-2 + H = hm 6 n = &oo- (31) 

The following theorem is a special case of a result stated by Herschfeld jHJ Theorem III] . For 
t — 2, it was first proved by Vijayaraghavan (in Appendix I of [121 p. 348]), and independently 
by Herschfeld. (The proof for t > 1 is similar to that for t = 2, and thus is omitted by both 
authors.) 

Theorem 11 (Vijayaraghavan-Herschfeld) Fix t > 1, and suppose that a n = a n (t) > for 
n > 1. The sequence (6 n ) n >i defined by (I30j) converges if and only if 



n— >oo 



Hmsupa n '* < 00. 
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1 /t n 

Note that this inequality is equivalent to the sequence (a n ) n >i being bounded. 

We now give the promised connection between the generalized Somos constant at and Ra- 
manujan's infinite radical. The latter is defined by 



f 1 + 2\/l + 3^yTT^ = lim B n (t) = B^t), (32) 

where 



B n (t) = ^ 1 + 2\J 1 + 3y 1 + • • • + (n - l)y 1 + n<fl. 
Corollary 2 t > 1 . // 



a„ = a n (t) = 1* -2* -3* •••n* 

/or n = 1,2,..., i/iera 

lim aj/*" = o\ < oo, (33) 

n— >oo 

and i/ie infinite nested radicals boo(t) and B^t) in (|31j) and f!32l) converge and are equal. 

Proof. The equality (J3*3*j) follows immediately from Definition |21 Since o\ < oo, Theorem ITT1 
implies that the infinite nested radical 



ai + V«2 + v^s + ■ ■ ■ = V I' 1 + V 142 ■ 2<1 + V 7 !* 3 • 2' 2 • 3* 1 + • ■ ■ = Jim 6 n (t) = &«,(*) 

converges. Now notice that o n (£) = -B n (t) for n, = 1, 2, Therefore, the infinite nested radical 

-Boo(t) also converges, and &oo(t) = B^it). This proves the corollary. □ 

Remark 4 The only known value of B^t) is -Boo (2) = 3, that is, 



l + 2yi + 3Vl + --- = 3. 

This formula was discovered by Ramanujan [T3J Question 289 and Solution, p. 323], but his 
proof is incomplete. Vijayaraghavan [TJJ p. 348] and Herschfeld [HJ pp. 420-421] each completed 
Ramanujan's proof. The result also appeared as Putnam problem A6 in 1966 (see |TJ pp. 5 and 
52]). 



4. Calculation of 7(2;) at Roots of Unity 

The purpose of this section is to calculate the value of 7(2) at roots of unity. We use the 
following summation formula, which has some resemblance to Problem 24 in |2Til p. 71]. 
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Theorem 12 Let (9 n ) n >i be a periodic sequence in C, and let q be a positive multiple of its 
period. If Ylj=i @j = 0> then 

n+1 A„, r ft 



V0 n log^ = V^log 



n=l j=l 

Proof. The Weierstrass product for the gamma function is [3 Section 8.322] 



00 z/n 



z AJ - 1 + i 

n=l n 

for z6C - {0, —1, —2, . . .}. It follows that for j = 1, . . . , g, 

logr f i) - io g r (ill) = 2 + log^tl + f; f--L + iog^±l±i) . (34) 

\qj \ Q J Q J q n qn + j ) 

Now multiply by 9j, and sum from j — 1 to g. Using X^=i = anc ^ = @qn+j, we obtain 

E f i log = E ^ log — + E E ° s log q -^ ±I±1 = E ^ ^g — , 

where the last series converges by Diri chief's test. This proves the theorem. □ 

Example 7 Take 9 n = (— l) n_1 and q = 2, and exponentiate the series. Using the identities 
(valid for z G C - {0, -1, -2, . . .}) 

T(z + l)=zT(z) (35) 

and 

r(z)r(i - z) = ^— (36) 

sin ix z 

to compute r(l/2)T(3/2) = tt/2, we recover Wallis's product for pi, which can be written 

An 2 7r 



n 



An 2 -l 2 ' 

n=l 

Now take 9 n = i n ~ x and q = A. Exponentiating the real and imaginary parts of the series, 
and putting z — 1/4 in we get the pair of products 

n f ( ~ 1)B " = r(1/4)2 fr = r (V 4 ) 2 



n=l x ' v n=l 
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Multiplying them together, or dividing the second by the first, we obtain the formulas 

r(i/4) 4 



x '2n + l^ ( - 1)n_1 



n 

n=l 



2n - 1 



n 

n=l 



4n 2 



7T 

T 



The last may be a new product for pi. If we multiply it by Wallis's product, the factors with 
odd n cancel, and taking the square root gives 



n 



71 



A 16n 2 - 1 2y/2' 
These products are perhaps known. 

The next result gives two formulas for the value of 7(2) at any root of unity u ^ 1. 
Theorem 13 Let p and q be relatively prime positive integers. If cu = e l7Tp ^ q 7^ 1, then 



7(w) + 



logfl - uj) 



00 



q 2q p ( 1+1 

^^-^ = ^^-Mog Vi/ 
3=1 3=1 



r ( 

'.29 



where 



E 



log 



j,p,g 



log- 



29 



r i±i r (2±a 



2? 



z'/p is even, 
if p is odd. 



(37) 



Proof. By Definition ^ and formula (jHJ), we only need to prove that 

00 q 2q p (i±± 

- E 1o k^ = E = E >° 1 



n=l 



3=1 



3=1 



rl .4 



(3f 



If p is even, (^'■ 7_1 )j>i is a periodic sequence of period q with ^jUi^ -1 = 0- I n this case , the 
first equality in (|3*5jl follows immediately from Theorem [T21 

If p is odd, (a; J_1 )j>i is a periodic sequence of period 2q with Xljli k^' -1 = 0- Hence in 
Theorem IT^l we may replace q with 2q, and take 6*^ = oji~ l . Since a;- 7 " 1 " 9 = — cu- 7 for j — 1, 2, . . . , q, 
the first equality in (|3*8j) holds with Ej >p>q given by (|3*7|) . 

To show the second equality in (|38|). use Theorem IT21 and the fact that a; is a 2gth root of 
unity. □ 
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Example 8 Take p/q = 1/2 and 2/3. Using Theorem El and identities (J33|) and (jSSjl . we get 



7(t) 
7 ( e *rf/3) 



1 /•! 



dx dy = — — log — ~ — — — ~~~ 

o Jo V 1 - ixy)(- log xy) 4 7r^/2vr 



+ Hog 



_8a/tt_ 

r(i/4) s 



1 — x 



o Jo (1 - e 2 W3 X y)(_l gxy) 



7T 3 , 2vr .1^/3.9 vr 

ax aw = — = log TT + i — log 

4^/3 2 6 3r(2/3) 2 I 2 6 2vr 12 



The following corollary gives an additional formula for the value of j(z) at any point ui of 
the unit circle with argument vrp/g, where p and g are integers with p odd and q > 0. For such 
a point u>, the formula expresses j(oo) as a finite linear combination of powers of uj with real 
coefficients. 

Corollary 3 Let p and q be integers, with p odd and q positive. If uj = e mp ' 9 , then 



where Ej tPtQ is given by (|3*7)l . and 



m 2q 



Proof. By Theorem 1131 and formula (jUJ), we only need to show that 



CO „_1 



E 



n=l 



3=1 



Let si, S2, . . • be the partial sums of the series. Then, using u 2q = 1, 

2<? fc-1 - 

lim s k = lim s 2g fc = lim V^ -1 V] — . 

j=l m =0 y J 

Since p is odd, uji +q = —uji for j = 1, 2, . . . , q. Therefore, 

g 



1 - 

lim Sfc = > a/ -1 — > 
i=i 



1 



m=0 [m+4-] (m+l±Z 



2q J y ' 2g 

Using the formula [UJ Section 8.363, Equation 3], 19, Equation 44:5:8 

1 i/j(x) - t/j(y) 



E 

m=0 



+ x)(m + t/) 



which is valid for x, y G C — {0, —1, —2, . . .} with x ^ y, equation (j4*Uj) follows. 



(39) 



(40) 



□ 
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Remark 5 In ()39|) the quantity if>(j/q) can be calculated from Gauss's formula 
< j < q, together with the value if)(l) = —7. 



when 



Example 9 Setting p — 1 and q = 3 in Corollary El we multiply by 6 and obtain 

6 7 (e^ /3 ) 



,'2\ , fi\ cl r (1/3) r (2/3) 
''Uj-n6j +61og r(i/6)r( 5 /6) 



+e 



iiz/3 



if, 



+e 



2mt/3 



-7 -if) 



. i\ , r (1/2) r (5/6) 
"'3j +61os r (1 /3) 



6 log 



6r (2/3) 



11V3 — 3 log + i 



7T 



7T 



r (1/2) r (1/6) 
/- r(i/6) 2 r(i/3) 2 y^ ' 

3^ log ^ 



The next theorem gives a formula for the average value of the function 7(2) at the vertices 
of a regular polygon inscribed in the unit circle with the positive real axis a perpendicular 
bisector of one side. 



Theorem 14 Given q 6 Z + , let u — e m l q . Then 

-1 ^1 



1 9-1 

^ fc=0 



1 — X 



1 

2^ 



Jo (I + xiyi)(- log xy) 
q + l 



dx dy 



2g 



— lot 



r U) r (f 

r (0 r (^) 



We give two proofs. 



Proof 1. By the method of partial fractions, for ( e C - {u 1 , u 



-{2q-iy 



I- 1 1 

^ 1 -(u 

k=0 * 



(uJ 2k+l 1 + C q 

Letting ( = xy and using (fTTj). we obtain 

9-1 »i »i 

/ i (l + x*yi)(-logxy) 



Y J l{^)=q 



1 — X 



dx dy. 



(41) 



A:=0 

Making the change of variables X = x q , Y = y q gives 

I- 1 rl /•! /w\±-l 



$> ( ^+. )= / / ^l'-^'Tj: 1 ^-./*-. 



fc=0 



^0 



;i + xr)(-io g xr) 
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Now use the following evaluations Corollary 3.4], which hold for u, v > 0: 



o 



'% . dxdY = i 

„ (i + A'y)(-iog.Yy) 2 



/o io (i + xr)(-iogXF) dAdy - M -^ og r(|)r(H±i)- 

Proof 2. Using Corollary El and the equality Ej t2 k+i,q = Ej t \ A from (|3T|). we obtain 

k=0 k=0 j=l j=l k=0 

Since to = e ln ^ q , if j £ {1, . . . , g}, then cj- 7-1 ^ —1, and = 1 if and only if j = 1. It follows 
that Efc=!)(y _1 ) 2fc+1 = when j G {2, . . . , g}, hence 

fc=0 

The theorem now follows using (jSZJ), (GUI), and (}4~Tj) . □ 

Example 10 Take q = 3. Using (J2~3~|) and (}3l)|) . we get 

7 (e^ 3 ) + 7 (e-) + 7 (e 5 -/ 3 ) ^ff 7— , dxdy = ^=- log 

Jo Jo (1 + s 3 r)(-loga;y) y/3 2 



5. The Hyperfactorial K Function and the Derivative of 7(2) 

In this section we study the derivative of j(z) at roots of unity, using the following function. 

Definition 3 The Kinkelin-Bendersky [TU1I2] hyperfactorial K function is the real- valued func- 
tion K(x) defined for x > by the relation 



logK(x) = — - -log27T + y \ogT(y)dy. 



(42) 



(Bendersky [2] uses the notation T ly but most authors today use K to denote the function 
defined by fl4*2^) ; see, for example, p. 135].) 

The K function satisfies [2, p. 279] K(0) = K(l) = K(2) = 1 and, for x > 0, 

K(x + 1) =x x K(x). (43) 
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Thus, by induction onfiG Z + , we have K{n + 1) = l l 2 2 ■ ■ ■ n n . 

There is an analog for the K function 2, p. 281] of Gauss's multiplication formula for the 
gamma function. A special case is 

n—1 / . \ , " 2 -i 

3 \ A n 



3=1 



n 



1 n^-l 



Jl 12n g 12n 



for n > 2, where pp. 263-264] 



A = lim f ^ + ^ 2 = lim 1 2 l22 """ n " = 1.28242712 . . . (44) 



n— >QO 



X 2 ^12 e 4 ^2^12e4 

is the Glaisher-Kinkelin constant (see P- 135]). In particular, 



K I 2 = (45) 



The following summation formula is needed in the proof of Theorem [TBI 

Theorem 15 Let (9 n ) n >i be a periodic sequence in C, and let q be a positive multiple of its 
period. If 

<? 

J> = 0, (46) 



then 



Proof. Dirichlet's test implies that the series converges. The proof of the formula for its sum is 
similar to that of Theorem The new ingredient is the Barnes G function, which is defined 
by the infinite product 

(2-kYI 2 -^r(l + -) n 



for zgC, and which is related to the K function by K(x + l)G(x + 1) = T(x + l) x for x > 
(see Section 2.15], Section 6.441]). 

From (|4*3*|) and ()35j1 . we can write the relation as i^(x) = G(x + l)~ 1 r(x) :c , so that 
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for j = 1,2, ... ,q. Now take the logarithm, multiply by 9j, and sum from j — 1 to q. We 
compute the result in two parts. On the one hand, a calculation using ()47|) and ()46|) yields 

i G ( i- + 1 1 i 5 oo<j ✓ . i • i i \ 
> # ? glog — V = / _ / > #7 ^ gnlog : — . 



On the other hand, from (|34j) we have 



£ id log -XU = - ± »,i i 1 - + log i±l) - £ ± e, Li. + j .og ^±i±l) . 

Therefore, by addition 

Substituting 9j = qn+ j on the right, and using (j4T?j) again, we deduce the desired formula. □ 

Example 11 Take 9 n = (— l) n_1 and q = 2. Exponentiating the series, and using (J43|) and 
fl43j) . we obtain the infinite product 



n 



n=l 



2V6 eV7r 



A 6 



where A is the Glaisher-Kinkelin constant 
Remark 6 The last equality, in the form of the limits 

2N+1 ( 2N ( 1 \(-V n+1 n A e 

lim ff 1 + - = e • lim TT 1 + - J = - T7?r - 

n=l v 7 n=l v 7 



is due to Borwein and Dykshoorn jSJ. 

Our final theorem is the main result of the section. It gives the value of the derivative of 
j(z) at any point uo 7^ 1 of the unit circle with argument a rational multiple of tt. 

Theorem 16 Let p and g be relatively prime positive integers. If uo = e mp ' 9 7^ 1, then the 
derivative of the function j(z) at z = uj is 



1 , log(l-u,) , * r (i)^(^) 2? 



2.3 



Proof. Using Definition^ the Taylor series of the derivative of the product z r y(z) when \z\ < 1 



is 



7 (z) + zi{z) = J2 ^ (l - n log — ) 

n=l ^ n J 



(48) 



For any fixed z / 1 with \z\ = 1, the partial sums of the series Yl^Li are bounded. It 
follows, using Dirichlet's test, that the series in (}4*H|) converges for such z. Since the function 
7(z) + z^(\z) is continuous at such z, Abel's limit theorem implies that (}4*%j) also holds when 
\z\ = l^ z. 

In particular, we may take z = to = e lnp ^ q 7^ 1. Then the sequence (u; n_1 ) n >i is periodic, 
and 2q is a positive multiple of its period. Since 5^n=i ujn ~ 1 = 0> m Theorem^Jwe may replace 
q with 2q, and take 9 n = uj n ~ l . Using Theorem IT^l the result follows. □ 

Compare the following double integral formulas with those involving A in [7] and [To*] . 

Corollary 4 TTie /zrsi and second derivatives of the function 7(2) a£ 2 = — 1 are 

„ r 1 r 1 xy(i-x) , , , 2 u / 6 a 6 

7(-l)=/ / 7; ; L r K dxdy = \og 



./o 



'1 + xy) 2 {— \ogxy) 7r 3 / 2 e 



and 

,„ 1N f 1 Z" 1 2zV(l-z) , , , 2 10 / 3 A 24 7C(3) 

7 " -1 = / / * v - ; — -dxdy = log , 1a/ . ^V. 

' l ; 7 7 (l + x U ) 3 (-logxi/) y 6 vr 4 e 13 / 4 2tt 2 

Proof. In both cases the first equality follows using Theorem El 

To obtain the value of 7 '(— 1) use Theorem ITol together with formulas (|43|) and (|45|) . Alter- 
natively, use pHjl. Example [TTJ and the value 7(— 1) = log-. 

To evaluate 7"(— 1), differentiate the inversion formula in Theorem |H1 twice at z = — 1 + zi, 
let t tend to + , and re- arrange terms, to get 



7 "(-l) = 4 T '(-1) - 3 7 (-l) + \ - 7 - logTr - j> Q) + Q ) + • 

Using the values Sections 44:8 and 44:12] 



vr V2/ 2 



the desired formula follows. □ 
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